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D3-XJENERATE  UNSTEADY  COMPRESSIBLE  FLOWS 


ABSTRACT 

Characteristic  conditions  are  derived  for  non-linearized  three  dimen- 
sional compressible  flows  in  which  the  entropy  of  any  fluid  particle  re- 
mains constant. 

Irrotational  isentropic  flows  for  which  the  maps  x,  y,  z,  t 

► d(//dx,  9^/ay,  dtjf/dz,  b<//dt 

have  fewer  than  four  dimensions  have  been  studied  by  applying  a degenerate 
Legendre  transformation  to  the  equation  of  the  velocity  potential  function, 
gl.  The  class  of  such  flows  includes  steady  conical  flows  and  pseudo- 
stationary flows,  for  example.  All  results  presented  are  in  the  nature  of 
statements  about  geometrical  properties  of  the  map,  of  certain  characteristic 
hyper surfaces  in  space-time,  etc. 

A number  of  problans  which  could  be  expected  to  lead  to  pseudo-station- 
ary flows  are  outlined.  Most,  but  not  all,  are  well  known  problems  that  have 
been  solved  at  least  in  linearized  form. 


1.  Introduction* 

The  development  of  the  theory  of  compressible  perfect  fluids  has  been 
impeded  by  the  fact  that  the  basic  equations  are  non-linear.  Therefore  it 
has  been  profitable  to  consider  special  examples  to  discover  general 
principles  and  to  gain  insight  into  the  structure  of  flows.  Specializations 
are  usually  based  on  such  properties  as  irrotationality,  reduction  in  number 
of  dimensions,  radial,  cylindrical,  or  axial  symmetry,  etc. 

For  steady  irrotational  flow  it  has  been  found  instructive  to  study 
those  flows  for  which  the  mapping  (x,  y,  z)-*-(u,  v,  w)  * V0,  where  0 is 
the  velocity  potential,  takes  the  physical  space  onto  a one-or  two- 
dimensional  image  in  the  hodograph  space  [j>J  This  leads  to  consider- 
ation of  the  class  of  conical  flows  which  are  interesting  and  important 
not  only  in  their  own  right,  but  alpo  for  the  fact  that  they  appear  to 
be  connected  with  the  singularities  of  steady  flow  fields  near  the  tips 
or  edges  of  bodies  with  discontinuous  normals*  That  is,  the  limits 
attained  by  the  velocity  components  and  other  flow  functions  at  the  tip 
will  depend,  in  genera}.,  on  the  direction  from  which  the  tip  is  approached. 
An  argument  can  be  devised  which  makes  it  plausible  to  conclude  that  the 
(u,  v,  w)  image  or  hodograph  is  the  same  as  that  for  the  flow  which  would 
occur  about  the  tangent  cone  at  the  tip  under  the  same  free  stream  con- 
ditions o 

Such  considerations  have  suggested  that  it  would  be  equally  profitable 
to  undertake  the  study  of  non-steady  irrotational  flows  for  which  the 
mapping  (t,  x,  y,  z)  -*~(30/at,  u,  v,  w)  » (00/dt ,V0)>  where  0 is  the 
velocity  potential,  takes  space-time  onto  a one-,  two-,  or  three-dimensional 
image.  The  image  will  be  called  the  hodograph,  and  when  it  has  fewer  than 
four  dimensions,  the  hodograph  and  the  flow  will  be  said  to  be  degenerate. 
The  study  leads,  eventually,  to  the  consideration  of  non-steady  conical 
flows,  also  referred  to  as  pseudo-stationary  or  self -similar  flows.  Such 
flows  occur  in  the  establishment  of  steady  conical  flows  about  conical 
bodies  in  shock  tubes,  not  to  mention  other  examples  discussed  at  the 
end  of  the  paper0  It  is  to  be  expected  that  they  are  also  related  to  the 
behavior  of  flow  fields  near  singular  points  in  space-time.  That  is,  the 
limits  attained  by  velocity  components  and  other  flow  functions  at  the 
singular  point  should  depend,  in  general,  on  the  direction  from  which  the 
singular  point  is  approached.  If  this  is  in  a region  of  potential  flow, 
one  would  conjecture  that  the  (80/st,  u,  v,  w)  image  or  hodograph  is  the 
same  as  that  for  one  of  the  degenerate  unsteady  flows  considered  here. 

This  phenomenon  would  occur,  for  example,  in  the  one-dimensional  flow 
behind  a piston  which  is  withdrawn  with  increasing  speed  and  which 
experiences  a discontinuous  increase  in  speed  at  some  instanto  The 
singular  point  would  be  on  the  piston  curve  at  the  time  and  place  where 
the  discontinuity  in  velocity  occurs. 

i.  Bracketed  numbers  designate  .references  at  the  end  of  the  paper0 
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Much  of  the  discussion  is  concerned  with  characht eristic  subspaces 
of  space-time  or  of  the  hodographs.  To  clarify  this  discussion,  and 
also  because  it  casts  some  light  on  the  general  wave  nature  of  unsteady 
flows,  a section  on  the  characteristic  conditions  for  general  unstea<$r 
flows  has  been  included, 

2,  Fundamental  Equations, 

Non-steady,  non-vi scous  flow  is  governed  by  the  equations  of  motion 

(2.1)  3u.  /3t  + u.3u ./dn./dx.  = - p "^Sp/ax. 

X 31  1 J X 

the  equation  of  continuity 

(2.2)  dp  /at  + a (p  u.)/ax.  = 0 

J J 

a pressure  density  relation 

(2.3)  p = es/cv  pY 

and  constancy  of  entropy  for  a fluid  particle  (between  passages  through 
successive  shock  fronts) 

(2th)  ds/dt  + u..as/ax.  = 0 

0 J 

Here  u^(x,  t)  (i  = 1,  2,  3)  are  cartesian  velocity  components  at  time 

t at  the  point  with  rectangular  coordinates  x^,  x„,  Xy  p(x,  t),  p(x,  t), 

and  s(x,  t)  are  the  pressure,  density,  and  entropy  per  unit  mass, 
respectively;  and  y is  the  ratio  of  the  specific  heat  at  constant  pressure 
to  that  at  constant  volume.  The  convention  has  been  adopted  that  repeated 
indices  in  any  sum  and  imply  summation  over  their  range.  The  pressure  can 
be  eliminated  from  (2,1)  by  means  of  (2.3).  It  will  be  assumed  hereafter 
that  this  has  been  done. 

Occasionally  it  will  be  more  advantageous  to  use  vector  notation. 

Then  the  preceding  equations  take  the  form 

(2.1v)  Su/dt  + u*Vu  = - p”1Vp 

(2.2v)  dp/dt  + V • (pu)  = 0 

(2.1iv)  as/at  + u*  Vs  =0 

In  some  parts  of  the  following  discussion  it  will  be  convenient  to 
consider  irrotational  isentropic  flows.  For  these  s **  constant,  and 

(2.?)  3 u.  /ax.  » 3u./3x. 

X J J X 


\ 
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Now  there  exists  a velocity  potential  function  <j>  such  that 

(2.6)  u.  = a ^/axi 

Then  (2.1),  (2,3),  and  (2.6)  imply  Bernoulli's  equation  dp/St  + t?  u^u^  + y p/(y-l)p 

df/dt  for  some  function  f(t).  Since  f(t)  can  be  absorbed  into  <$  without  affecting 

(2.1)  to  (2.6),  Bernoulli's  equation  can  be  expressed  as 

(2.7)  a2=  9p/Sp  = y p/p  = - (y  - l)(-jj  u^u^  + a0f/at) 

where  a is  the  speed  of  sound.  Finally,  eliminate  p and  p from  (2.1)  and  (2.2) 
by  meaEs  of  (2.3)  and  (2.7)  to  obtain 

(2.8)  (a2 6.  . - u.u.)  3 (2f/3x. 3x . - 2u.a?5f/ax  at  - a2(2f/t2  * 0 

-L  J J J J 

where  Kronecker's  delta,  6. . = 1 (0)  accordingly  as  i = (/)  j. 

J 

3.  Characteristic  Conditions. 

A standard  exercise  in  the  theory  of  partial  differential  equations  is 
Cauchy's  problem.  For  the  system  (2.1v),  (2.2v),  and  (2.hv)  this  can  be 
stated  as  follows.  Consider  a three-dimensional  initial  hypersurface 
in  soac e-time. 

(3.1)  x = T(r1#  r2,  r3)  = X(r),  t = T(rx,  rg,  r3)  » T(r) 
where  the 

(3.2)  rank  of  \\  aY/ar^  dH/dr±\\  = 3 

Construct  solutions  u(x,  t),  p(x,  t),  and  s(x,  t)  which  assume  on  (3.1)  the 
prescribed  initial  values 


u(x(r ), 

T(r ) ) = TJ(r) 

(3.3) 

p(I(r), 

T(r))  = R(r) 

s(Y(r), 

T(r ))  = S(r) 

subject  to  appropriate  assumptions  about  the  continuity  and  differentiability 
of  the  functions  (3.1)  and  (3.3). 

If  the  function  T(r)  is  a constant,  TQ,  then  (3.1)  is  merely  some  three- 
dimensional  region  of  the  x^x2x3  space  at  a fixed  time  t **  TQ. 

If  T(r)  is  not  a constant  function,  then  for  the  solutions  of  T(r)  a TQ 
= constant  the  functions  X(r)  define  a surface  in  the  x-jX^- space, 


7 


the  shape,  location,  and  orientation  of  which  depend  on  the  time  TQ. 

In  other  words,  (3.1)  can  be  interpreted  as  defining  a moving  twisting, 
deforming  surface  in  the  x^XgX^-space. 

If  1,  T,  If,  R,  and  S are  analytic  in  r^,  r2,  and  r^,  try  to  expand 
u(x,  t),  p(x,  t),  and  s(x,  t)  in  Taylor's  series  about  some  point  I*,  T* 
of  (3.1).  The  values  of 

(3.4)  Vu,  au/at,Vp,  ap/at,  Vs,  as/at 

required  for  this  expansion  must  satisfy  (2.1v),  (2.2v),  (2.4v),  and 
the  strip  conditionsi 

(3.5)  \ *Vu  ♦ T^/at  * \ 

(3.6)  \ °Vp  + ^ap/at  = 

(3.7)  \ *ys  + T^s/at  - 

where  fi  * df/br±  for  any  f(r),  and  i ■ 1,  2,  3.  If  these  twenty 

equations  have  a unique  solution  (3.1*)  on  (3.1),  then  the  values  of 
all  partial  derivatives  of  all  orders  are  uniquely  determined  on  (3.1). 
By  the  cauchy-Kowalewski  theorem,  the  Taylor's  series  for  H(X,  t), 
p (3E,  t) , and  S(x,  t)  are  unique,  converge  near  X#,  T*,  and  satisfy 

(2.1v) , (2.2v),  (2.4v),  and  (3.3)‘. 

In  all  flow  problems  u,  p,  and  s must  be  continuous  in  regions  of 
space-time  between  shock  fronts  or  contact  discontinuities,  but  (3.4) 
or  the  partial  derivatives  of  higher  order  merely  need  to  be  continuous 
almost  everywhere,  so  tr,  p and  s are  not  necessarily  analytic.  Hence 
some  partial  derivatives  of  u,  p,  or  s may  have  a discontinuity.  If 
this  is  interpreted  as  a weak  disturbance,  it  is  important  from  physical 
and  computational  considerations  to  determine  how  this  is  propagated. 

For  simplicity,  suppose  some  of  (3.4)  are  discontinuous  at  J*,  T*  on 

(3.1) ,  though  the  discussion  can  also  be  adapted  to  partial  derivatives 
of  higher  order.  Then  it  must  be  impossible  to  determine  (3.4)  uniquely 
at  Y#,  T*  from  (2.1v),  (2.2v),  (2.4v),' and  (3.5)  to  (3.7),  but  this 
system  must  still  be  compatible.  For  theoretical  and  numerical  appli- 
cations the  most  interesting  situation  is  that  in  which  (3.4)  fail  to 
be  uniquely  determined  everywhere  on  (3.1)  rather  than  just  at  I#,  T*. 

In  this  case  (3.1)  will  be  said  to  be  a characteristic  hypersurface, 
and  the  compatibility  conditions  for  the  system  of  twenty  linear 
equations,  which  are  expressed  solely  in  terms  of  (hyptr-) surface  data 

(3.1)  and  (3.3)  and  their  surfacq  (tangential)  derivatives,  will  be 
called  characteristic  conditions. 


An  application  of  Cramer »s  rule  to  (3.1*)  yields  (Xj.  X^Vu  - 

^j(Xj+iXXj+2)(^  - Tj3u/dt)  s where  (X  1 U)  *X«1?xS  for  any 

vectors  X,  Bs  C,  and  j + 1 and  j + 2 are  to  be  reduced  mod  3.  Now 
the  four  components  NqS  N of  a normal  to  (3.1)  satisfy 

(3.8)  NQT±  + JMj.  - 0 

By  (3.2)  this  system  has  one  linearly  independent  solution  :NQ  ■ (X,  i 


s »- 

Hence 

N V u 
0 

=*  lau/at 

(3.9) 

»0Vp 

- Iap/at 

noV. 

= Nas/at 

* 

Then  by  (2.1v)s 

(2.2v)s  ( 

2.1*v)s  and  (3.9) 

(3.10) 

(»„* 

TMOdS/at 

+ R"1  N [k2dp/dt 

♦ (ap/a  s)as/at]  + 

2 

3<V* 

X>2>  L*™j  i »'  < 

(a2r_p  s^ap/as)]  - o 

(3.11) 

(N0« 

ii*N)ap/a 

t +■  RN*au  t + Z.(xi+1x  x.+2  )*(rTJ)  - o 

(3.12)  (no  + u»N)as/a  t +2^  Xj+1Xj+2)Sj  - o 

where  P * eS/CvRY  and  A2  ■*  8P/BR.  The  scalar  product  of  and  (3.10)  yields 
(NQ  + ^l?0N*eiu/9t  + R~1N2  |^A29p/at  + (dP/as)ds/3t]  + 

V Sj^P/as) (V Iw  5^)]  . 0 

The  determinant  of  the  coeficients  of 
(3.11*)  N»3u/ats  ap/ats  as/at 

in  (3.11)  to  (3.13)  is 

(3.13)  A » (No  + fl-I)  [(N0  + flVft)2  - A2N2] 

9 
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If  A 5^0,  then  (3.11i)  are  uniquely  determined,  and  since  NQ  + TMJ  f 0, 
3u/3t  are  unique  by  (3*10)  If  in  addition  N^O,  then  by  (3.9)  Vu.,  Vp>, 
and  Vs  are  unique.  Thus  a necessary  condition  for  non-unique  values 

of  (3*10  is 

(3.16)  NqA-0 

Case  Is  If  the  Jacobian 

(3.17)  N0  * Xg  X3)  - 0 

its  rank  is  two  by  (3.2),  and  also  T^O  for  some  i.  Hence  F(5)  * 0 for 

some  F.  In  other  words,  the  initial  surface  in  the  - space  does 

not  vary  with  th|  passage  of  time.  New  (3*9)  are  compatible  since 
N*X^  » 0,  and  hence  N X * °»  Their  unique  solutions  satisfy 

N2^u/at  - 2 j(W  *j+1  Xj+2)TI.3  N2dp/3t  * ^(fl  IJ+1  xj+2)Rd 


»2a.M  - 2.  IS  J ^ 1^)3 j 


By  means  of  these  equations  eliminate  ciu/dt,  ap/at,  and  Ss/at  from 
(3.10)  to  (3«12)  tq  obtain 

2 j(Tf  ^+1  Xj+2)  [tMI^  ♦ + SjdP/dB)  ] + 


(3.18) 

(3.19) 

(3.20) 


IT2  ♦ r“1(a2r^  + s^ap/as)  ] - o 

^>1  *J+2^N  + " 0 


■3 


Xj+1 1^2)  + n20J  Xj+1  I>2  )]s3  - o 


Case  2s  Suppose  NQ^0,  but 


(3.21) 


n + U*H  * o 
0 


Then  the  compatibility  conditions  for  (3.10)  become 

(3.22)  ^X^^j+lXXj+2^  [>,  +H"1-(A2R;j  ♦ Sj3P/aS)]  - 0 


10 


If  these  are  satisfied,  then  (3.10)is  equivalent  to  the  single  equation 
(3»13)j  'which  now  takes  the  form 

(3.13.2)  0 =*  R_1N2  A23p/at  +(0P/9S)8s/at  + 

3 *3*1  V’  + R‘1(a2r3  * sJ8p/i>s,(*  V V> 

while  (3*11)  becomes 

(3.11.2)  RU»auM+  3 0 

Finally  (3.12)  yields 
(3.23)  “j+2)Sj=  0 

By  (3.11.2)  and  (3.13.2)  three  of  the  five  functions  (3.lU)  are 
arbitrary.  Hence  there  is  a three-parameter  family  of  solutions  (3.U). 

To  interpret  (3.21)  proceed  as  follows.  Suppose  a solution 

u(x,  t),  p(x,  t),  s(x,  t)  of  (2.1v),  (2.2v),  and  (2.Uv)  is  known.  First 
seek  an  integral  t = f(x)  of  (3.21).  In  this  case  N = - 1,  N * Vf  = p, 
and  (3.21  becomes  ' ° 

(3.2U)  u-p  - 1 =»  0 

where  u = u(~,  f).  Solutions  of  this  first  order  partial  differential  1 
equation  can.be  constructed  by  passing  through  each  point  of  some  two- 
dimensional  initial  manifold 

(3.25)  x = x*(r1,  r2),  t = t*^,  r?) 

integral  curves  of  the  system 

cSx/dr ^ = u,  1 dt/dr^  =»  u*p  = 1 

(3.26)  dp/dr^  * - ( u *p  - (p’Bu/atJp 

with  initial  conditions  u = u*(r^,  r2)(known,  of  course)  and"p  =»  p#(r^,  r2) 

chosen  to  be  solutions  of  u#  *p»  = 1 and  3t*/9r  + p*  *ax*/3r  = 0 

ct  a 

(a  = 1,  2),  By  (3.26)  assume  r-  = t with  no  loss  of  generality.  Hence 
dx/dt  » u.  In  other  words,  at  all  times  t the  surface  T “ x(rp  r25 

is  always  composed  of  the  same  particles,  i.e»,  the  integral  surfaces 
t«f(x)  of  (3.21)  are  particle  fronts . 


m 


♦ 


\ 


n 


31 

Integral  surfaces  of  (3.21)  which  cannot  be  represented  by  t = f(x) 
must  satisfy  some  equation  g(Tc)  » 0,  where  g is  independent  of  t.  On 
such  a siirface  * 0,  contrary  to  the  hypotheses  for  this  case. 


Case  3:  "Suppose  NQ(NQ  + TTTJ)/  0,  but 


(3.27)  (Nq  + TH7)2  - A2N2  = 0 


By  (3*12)  as/at  is  unique, 
(3.13)  is 


A compatibility  condition 

(NQ  + 'U*TT)T!-I^  - R_1N2(aV. 


for  (3.11)  and 
+ S,9  P/9S) 

J — 


+ 


(3.28) 


anrj+1  Xj+2)  (N0  + U-Dr^CaV  + s^a  ?fas)  - 


a2:t2(TJ . x.  . I.  0 ) 

0 3+1  3+2 


}■ 


If  this  is  satisfied,  only  three  of  (3.10)  and  (3.11)  are  independent,  and 
there  is  a one -parameter  family  of  solutions  (3.U). 

To  interpret  (3.27),  again  suonose  a solution  of  ufx,  t^,  p(x,  t), 
s(x,  t)  of  (2.1v),  (2.2v),  and  (2.hv)  is  known.  Seek  integrals  t «*  f (x) 
of  (3.27).  Again  N = - 1,  Jf  = V f -o,  and  (3.27)  becomes 
(3.29)  (u-p  - 1)2°  - a-  p2=  0 


Solutions  of  this  first  order  partial  differential  equation  can  be 
constructed  by  passing  through  each  point  of  some  two-dimensional  initial 
manifold  (3.25)  integral  curves  of  the  system 


(3.30) 


dx/dr^  -(u-p  - lju  - a2p,  dt/dr^  = p*dx/dr^  = u-p  - 1 


dp/dr^  = - (u.p  - 1)  (Vu).p  + p2  a Va 
- [(u.p  - 1)  ut-p  - aafcp2 


with  initial  conditions  p*(r^,  r?)  determined  by 

(u*.p*  - l)2  - a*2  P*2  = 0 
3t*/3r  + p».9x*^r  =0,  a = 1,  2 


(3.31) 

(3.32) 
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An  important  special  case  is  that  for  which  (3.25)  degenerates  into  a 
single  point  x*  = constant,  t#  » constant.  Then  (3.32)  need  not  be 
considered,  and  (3.31)  yields  a two-parameter  set  of  initial  values 
p*  (r1,  r2).  The  set  of  all  integral  curves  of  (3.30)  through  x*,  t* 

for  each  set  of  initial  values  forms  a characteristic  hyperconoid. 

Eliminate  r^  from  (3.30)  to  find 

(3.33)  dx/dt  = u - a2p/(u-p  - 1) 

and  a similar  equation  for  d!p/dt.  At  each  instant  t,  the  points  xfr^rg  J t) 

on  the  characteristic  hyperconoid  integrals  of  (3.33)  determine  the_boundary 
of  the  region  into  which  a weak  disturbance  at  (or  sound  wave^from  x*  at 
time  t#  has  propagated.  By  (3.29)  and  (3.33)  (dx/dt  - u)  = a ♦ Thus  to  find 
the  possible  initial  velocities  dS/dt  at  time  t#,  to  draw  all  vectors  from 
the  point  x*  in  XjXgX^-space  to  the  surface  of  the  sphere  with  center  at 

x*  + u#  and  radius  a-::- . This  yields  the  well  known  result  that  weak  disturb- 
ances at  5c*  at  time  t*  propagate  in  all  directions  in  subsonic  flows,  but 
only  jtvEo*The  downstream  nappe  ofThe  Mach  conoid  with  vertex  at  x»  in  super- 
sonic flows. 

To  put  the  results  just  obtained  into  a more  familiar  light  consider 
one-dimensional  unsteady  flow.  Now  (3.1)  can  be  given  the  form  t «*  t(r^), 

x = (x(r^),  r2,  r^),  and  u = (u(r1),  0,  0).  (3.21)  and  (3.23)  yield 

(3.33)  dx  *’udt,  ds  = 0 

while  (3.27)  and  (3.28)  yield  (dx  - udt)2  = a2(dt)2,  and  if  dx  / 0, 

(dx  - udt)  du  + p dt(a  dp  + 9p/as  ds)  = 0;  In  their  usual  form  these  are 
(3.35)  dx  = + adt,  du  = + ^(adp  +a  "''9p/9sds) 

Similarly,  for 

2 2 2 

u = (ux/r,  uy/r,  0),  r = x + y 

2 2 2 

u = (ux/r,  uy/r,  uz/r),  r=x  + y + z J 
where  u ■ u(r,  t),  and  for  (3.1)  of  the  form  t ■ t(r^), 

x » (r(r^)  cos  r2,  r(r^)  sin  r2,  r^) 

3 (rC^)  cos  r2  cos  Tyt  r(r^)sin  r2  cos  r^,  r(r^)  sin  rj 

equations  (3.21),  (3.22),  (3.27),  and  (3.28)  reduces  to  the  dliaract eristic 
equations  for  unsteady  cylindrically  £ spherically  J symmetrical  flow. 


\ 
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U.  Specializations.  Degenerate  Legendre  transformations . (2.1),  (2.2), 
and  (2,4)  form  a system  of1  non-linear  partial  differential  equations  for 
five  functions  of  four  independent  variables.  At  present  the  theory 
is  much  more  highly  developed  for  systems  with  two  independent  variables 
than  for  three  or  four.  Hence  it  is  profitable  to  consider  examples  of  un- 
steady flows  with  special  properties  that  make  it  possible  to  reduce  the 
number  of  independent  variables,  A common  procedure  is  to  consider  one- 
dimensional or  cylindrically  or  spherically  symmetrical  fLows,  A slightly 
more  general  technique  which  has  proved  fruitful  in  the  study  of  steady 
flows  is  to  seek  solutions  of  the  form 

(4.1)  u±  * tu(|j.),  p - p (|X ),  s * s(|x) 

where  ji  =*  ••«,.  pnj  pa  = x2*  x3*  ^ = “ 1>  •••»&» 

1 n«  4j  and  the 

(4.2)  rank  of  ||  3pa/9xi,  3pfl/3t  ||  ■ n 

As  in  the  steady  case,  this  leads  to  very  interesting  classes  of  flows. 

Several  important  examples  of  this  type,  investigation  of  which  should 
yield  valuable  results  both  in  fluid  mechanics  and  in  the  theory  of  partial 
differential  equations,  will  be  mentioned  in  Section  8.  For  simplicity 
the  discussion  in  this  and  the  three  following  sections  will  be  restricted 
to  irrotational  isentropic  flow.  Now  s = constant,  and 

(4.3)  ut  = 3^/fexi 

where  p,  and  p satisfy  the  partial  differential  equation  (2.8)  and 
Bernoulli's  equation  (2.7).  Let 

(4.4)  xq  * t,  uQ  =*  3^/3t 

where  by  (2.7)  uQ  = uQ(p).  Assume  that  the 

(4.5)  rank  of  ||  3uI/3pa  ||  - n, 

where  I * 0,  1,  2,  3. 

By  analogy  with  the  usage  in  , the  image  u^  of  Xj.  under  the  mapping 
Xj.— > Uj.  will  be  called  its  hodograph.  Now  (4.1)  defines  an  n-dimensional  sub- 
space in  the  hodograph  space.  The  flow  will  be  said  to  be  a single,  double, 
or  triple  wave  accordingly  as/ft-i,  2,  or  3J  n-dimensional  subspaces  oT 
either  the  hodograph  space  or  of  space  time  will  be  called  curves  for  n = lj 
surfaces  for  n = 2 and  hyper  surfaces  for  n * 3.  If  n ^ 3 the  hodograph  and 
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the  flow  will  be  said  to  be  degenerate,  With  no  loss  of  generality  assume 
that  1=  n s!3.  Now  let  k(x,  t)  « # - Xj-u^p).  Then  ak/fexj  = - 

(Xj.3Uj.^p,o)  (3^/aXj),  and  the  rank  of  ||  ak^Xj,  3pa/cbCj  || . * n also,  which 
implies,  k = k(p),  Now 
(tu6)  4 " k(p>  ♦ XjUj^) 
and  by  ( k,2 ) 

(U*7 ) a^BUj/ap,  + ak/9(ia  =*  0 
whence 

(U.8)  (xj  a^/a^a^p  + a^/apgap^app/kxj  * - a Uj/apa 
Again  by  (U.2)  the 

(U.9)  rank  of  |j  Xj  a2  u^/ap^pp  + B2k/3pa  dppfl3  n 

Hereafter,  assume  that  Uj.(p)  and  k(p)  have  been  chosen  to  satisfy  (lu9). 

Than  (U,7)  can 'be  inverted  to  yield  pQ(x), 

By  (Ini)  a point  on  the  hodograph  is  determined  by  setting  pQ  * PaQ, 

The  set  of  points  in  space-time  which  is  mapped  onto  Uj(pQ)  will  be  called 
its  prototype.  .(In 7)  implies 

Theorem  U.1:  Suppose  the  coordinate  axes  of  space-time  and  the  hodograph 

spaces  are  parallel.  Then  the  prototype  oT" a point  ^(p0)  on  an  n- dimen  - 

sional  hodograph  is  the  (U  - n)  - dimensional  intersection  of  n hyperplanes 
(In 7)  with  normals  parallel  to  the  tangent  vectors  3 Uj/9pg  at  tij.  (pQ;, 

Let  Dq  = (aui/3p(j ) (aui/dpQ ) , where  a is  not  summed.  If  = 0,  but 
d uo/app  i 0 for  seme  p,  then  (in  7)  defines  a (1+  - n)  - dimensional  region 
of  XjXjX^- space  at  time  t a 1 (ak/Sp^ )/ (Bu^p^ ),  If  0 for  all  c , then 
(b*7)  corresponds  to  the  intersection  of  a set  of  n planes  in  x^xy- space 
which  move  with  constant  velocities  - (9u  /fep  )(3u.  /Sp  )/D  , where  a is  not 
summed.  That  is  the  points  x^  that  bear  given  values  u^fu.^).  p(pQ),  and 
p Cli0>  move  a constant  velocity  that  depends  on  paQ,  This  suggests 


\ 


1* 


examining  the  possibility,  considered  at  length  in  Section  8,  of  seeking 
flows  with  uniformly  expanding  geometry,  in  which  all  flow  functions 
depend  only  on  x. /t . so-called  unsteady  conical  flows  [2J  or  pseudo- 
stationary flows1  fh,6l  . - " 

So  far  $ has  merely  been  compelled  to  yield  a degenerate  map.  For 
an  unsteady  compressible  flow  (2.7)  and  (2.8)  must  also  be  satisfied.  In 

(2.8)  32  (2f/9Xj.9Xj  is  required.  By  (iu3)  and  (U.I4.) 

(U.10)  bV/^Xj-SXj  - (auj/0^)  Opa/axj)  - (3Uj/a,xa)  (a^/fe  Xj) 

2 

Let  Vjp(p,)  be  U - n mutually  orthogonal  unit  normals  to  Uj.  ■ u^  (jx) • 
That  is, 

(U.ll)  = 0,  vIfvlA-  6rA 

By  (lulO)  (aUj/a  jaq)  (3pa/3  Xj)vj^=  0,  and  by  (lu5) 

(U.12)  (Sp^/c)  0 

By  (lu5)  and  (lull) 

(li.13)  3na/0Xj  = AapauJ/9^ 

for  some  A^.  Now  by  (H.8),  (lulO),  and  (iu£)  ( x^dTu^/d 2Uj + 

3 k/3p.a3pp  ) * - 6^.  The  general  solution  of  (h.7)  is  of  the  form 

(li.lli)  Xj.  - ^ + rrvir 

where  x^Q  is  a particular  solution,  and  rj_are  independent  of  the  p,'  s 'and 
each  other.  Write 

(k,l$)  - Afl  3 Uj/3^ 

2.  To  summarize  the  conventions  regarding  indices  in  the  following 
discussion,  Roman  lower  case  range  from  1 to  3;  Roman  capitals  from 
0 to  3;  Greek  lower  case  from  1 to  nj  and  Greek  capitals  from  1 to  k - n. 


5®#isS 
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Then  by  (U.7) 

(1.16)  Aa^3ul/9^o^  (an,-/9|J.p)  + 9k/9|ip  - 0 

defines  A (p),  and 

(1^.17)  [(^  + rr  vir  + 82k/3  i^Spp]  A^  - - 6ay 

defines  A^  (tx,  r).  Plnally,  by  (2.8),  (U.10),  and  (U.13) 

{f  (3  - (3  uq/3  p,a)  (3  %/fyip)]  - 

(U.18) 

0 V9^a  + ui9  (S  V%  + V V^p^ap  “ ° 

5.  Simple  waves . When  n * 1,  AQp  consists  of  a single  element,  and 
(U.18)  yields 

(5.1)  (u^  + UjUJ)2  = a2ulur 

where  primes  denote  ordinary  derivatives  with  respect  to  and  by 

Bernoulli's  equation 

2 X 

(5.2)  a - - (y  - 1)  (uQ  + ) 

If  a2  4 0 (5.1)  and  (5.2)  imply 

(5.3)  u!u|  » Ua,2/(y  - l)2 
Now  let 

1 

(5.U)  dtf/d^  = (u*  u,)?  =*  + 2a' /(y  - 1) 
to  obtain 

(5.5)  o + 2a/(v'-  1)  = o0  + 2a0/(y  - 1) 

where  aQ  and  aQ  are  constants.  Note  that  a is  simply  the  arc  length  of 
the  projection  of  Uj.  « ^(p)  onto  any  hyperplane  uQ  ■ constant.  In  a 

one-dimensional  simple  wave  (5.5)  reduces  to  the  familiar  form  (3*25) 
with  s ■ constant.  Since  (5.1)  is  of  the  form  (3.27),  the  prototype 
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hyperplanes  are  characteristic  hyper  surf  aces.  Since  u. , p,  and  s are 
constant  thereon,  (3.28;  is  satisfied,  further  note  that  each  proto- 
type hyperplane  corresponds  to  a plane  moying  with  const ant^velocity 

“ 7 + a]  uj  (u*u' )“  . 

— J is  m nr 

_ 1 

Since  u^ul(u^u^)-  "2  is  the  component  of  fluid  velocity  normal  to  this 
plane,  its  speed  is  exactly  what  one  would  expect  in  one-dimensional  flow. 

In  general  k(cr)  and  two  of  the  functions  u^cr)  can  be  chosen  arbi- 
trarily in  constructing  a simple  wave  in  accordance  with  (U.7),  (5.2), 
and  (5.5). 


6.  Double  waves.  When  n ■ 2,  let  g 0 be  the 

-I.  .1 ...  — ■ — . 

the  surface  Uj.  * ^(p.)  based  on  the  euclidean 
\b  be  two  second  fundamental  tensors  . 


covariant  metric  tensor  of 
2 

metric  ds  ■ du^du^,  and  let 
By  definition 


(6.1)  gap  « (9  Uj./9(io)(9  Uj/b  p,p) 

(6.2)  bfp  * vM€  a2  ityfe  p.a  9 (Up 

where  are  two  mutually  orthogonal  unit  normals  which  satisfy  (lull). 


Also 

(6.3)  Uj^./a  (xa  «*  qaq/a  p.a 


p 

where  <1  = \ Write  9q/9(j.a  «*  q,a 
covariant  derivatives  with  respect  to 


where  the  subscript 

,a 

|xa  and  based  on  gffp 


denotes  the 


By  (U.15),  (lul 6),  and  (6.1) 


(6.it)  * - (9  k/kjj^  )gr6(9  UjAig) 

where  gY6  is  the  inverse  of  g^.  Since  the  Christoffel  symbols  of  the 
first  kind,  based  on  gap,  are  [a  p,  yj  * (9  ^/9|iy)(92^/^a9^p%  the 
second  covariant  derivative  of  k becomes 


(6.5)  k^p  * 92k/9(ia9|ip  - (9k^nY)gY5(9 

Now  (h*  17) becomes  (k  flp  + r£bap)Ap^“  " ®ay*  S^nce  (^*18)  is  homogeneous 
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in  Aap,  only  the  adjoint  of  (k^  + - (-DC+P(^a+1  p+i  Cl  p+l> 

is  required*  Then  since  all  four  of  r6  and  ^ are  independent,  (U.18) 
yields 

(6.6)  [a^ ( gc p ”uo , auo , ]3 ) “ K,a+*V)(  (-1)a+P  Cl  p*l  “ 0 

.7)  [i 


(6 


a .«)  “ K,a+*V(%,P^V 


'ap  o,a  o,p 
Note  that  the 


(-  if+\ 


ia+1  P+l  “ 0 


(6.8)  rank  of  ||  a (ga?  - u^u^)  - (un  „ +qq  „)(un  R + qq  fl)  |j>0 


ft  n O T HU  ft, 

o,a  ,a  o,p  ,p 


For,  suppose  the  rank  were  zero.  Then  uq  Q + qq  Q/  0 for  some  a since  gQp 

is  non-singular.  Let  L be  a non-trivial  solution  of  (u  + u.u.  )L.  « 0. 
^ a o,a  x x,a  a 

Now  rank  zero  implies  (u.  L )(u.  aLQ)  3 0,  whence  u.  L 3 0.  Thus  u. 

x,<xe  x,[3p  x,a  a x,o 

A.  (u  + u.u.  ) for  some  A,-  * If  u.A.  3 0,  then  u.  =»  A.u  . and  g „ ■ 
x o,a  j jjU  x x x x,a  a.  o,a'  ap 


(1  + A-A^)  uQ  auQ  p,  which  is  singular.  If  u^A^  * 1,  then  uQ  a - 0,  so  UQ  « 
30/3t  = constant,  and  ^ = uQt  + f(x^,  Xg,  x^)  which  yields  a steady  flow  field, 
to  be  excluded  from  the  present  discussion.  If  u.A.  / O or.l,.  u. 

■ Ai  %c/(1  - ttjV«  Wt'en0e  (1  • Vi>  g«P  "L(1  - Vi>  * AiAl]uo,auo,fi 

is  singular. 


(4,4)  is  a system  of  two  second  order  quasilinear  partial  differential 
equations  for  four  functions.  To  determine  Uj.(p)  requires  two  more 
equations  which  may  be  obtained  by  assigning  "a  special  form  to  the  co- 
efficient tensor  of  (6.6).  The  resulting  systems  are  classified  according 
to  the  nature  of  the  characteristic  curves  of  their  integral  surfaces,  A 
characteristic  is  a curve  "on  Whxdh  the  coordinate  functions,  their  first 
partial  derivatives,  and  hence  the  metric  tensor  are  continuous,  while  the 
components  of  the  second  fundamental  tensors  may  have  discontinuities. 

Suppose  Bu^/Sp^  are  known  along  pa  = jxr(t)  on  u^  » Uj.(p),  By  (6.2)  the  strip 

conditions  d(  3x^/3^  )/dt  ^9%/%3pp)d|i£/dt  imply 


(6*9)  ‘ bap  * vM€d^a 


Then  b^  fail  to  be  uniquely  determined  along  pa 
if  and  only  if 

(6.10)  [a2(g(tf  - Vo),j)  - (u0jI*  l<l)a)(u0jp* 


* |ia(t)  by  (6.6)  and  (6.9) 
qqjP)]  (dp.a/dt)(dtip/db)  - 0 


\ 
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This  defines  the  characteristic  directions  d^/dt.  By  (6.5),  if  Uj(n)  are 

known,  then  (6,7)  is  a linear  partial  differential  equation  for  k which 
also  has  the  characteristic  directions  (6,10).  Equations  (6.6)  and  (6.7) 
will  be  said  to  be  of  hyperbolic,  parabolic,  or  elliptic  type  whenever  fl  * 

drt  ||a2(gop  - u0(<[  u0>p)  - (u^  * M)0)(u0;p  * <N>p)||<0>  ‘ °>>0-  Herc- 

after  (6.6)  and  (6.7)  will  be  assumed  to  be  hyperbolic. 


If  ac  is  the  arc-length  of  a characteristic,  then  (6,10)  becomes 

(6.11)  (du0/cbc  + qdq/d<rc)2  = a2  [l  - d^/cbj2] 

which  is  identical  with  (5.1).  Hence 


"heorem  6.1s  The  characteristics  of  the  hodographs  of  double  waves  are  one- 
dimensional holographs. 


By  (6.11)  the  component  of  u + qq  along  either  characteristic 

r ia  1 

is  (uo,«  + qq,a}  <VdcTc  “ 1 a L1  **  (duo/dorc)  J?  * Hence 


Theorem  6.2:  On  the  hodographs  of  double  waves  the  curves  of  constant 
speed  of  souncTand  feeir  orEhogon,- 1 trajectories'~Ersect  the  angles  between 
the  character!  stTcs. 


In  general,  curves  x^.  = Xj.(t)  (where  t is  any  parameter,  not  necess- 
arily time)  in  space-time  are  mapped  onto  curves  jx  = pa(t)  on  the  hodograph. 

It  is  convenient  to  know  the  relation  between  tangent  vectors  of  a pair  of 
corresponding  curves.  By  (U.7),  (lull*),  (6.U),  and  (6.5)  (6,12) 

(6.12)  (dxj./dt)^/^)  + [r£(t)ba^  + kjaJ  d^/dt  - 0 

for  some  r,  ( t ) . Unless  det  ||  r€  b^  + k j = 0,  this  determine-  dp^/dt. 


Conversely,  if  the  curve  p.^  » p.^  ( t ) is  given  on  a hodograph  surface,  its 
prototype  is  the  hypersurface 

(6.13)  ^(r,  t)  * x^t)  + jr€  - A£(t)J  vIg(t) 

where  the  normals  v^to  the  hodograph  satisfy  (U.ll),  Xj.Q(t)  is  defined 
by  (6.U),  and  A^(t)  by  dA£/dt  * v.^  dx^/dt  (to  ma!<e  the  three  families  of 
c'irves  of  naratieters  r^,  or  t mutually  orthogonal,  as  a matter  of 
convenience);  For  (6,13)  an  analog  of  (6.12)  is 

(6.1U)  (axj./BtJOuj./a^)  + [(r€  - k€  )bap  + k dp.p/dt  = 0 
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Since  vTgdXj./dt  * 0,  (6.1U)  implies 

(6.15)  axj.  /dt  « - (3uI/Dixft)gap  (r6  - A4  )ba^  + k^  dj^/dt 

In  general  the  direction  of  Sx^/Bt  varies  with  r^  and  rg  on  the 

prototype  of  p,  (t),  so  the  hypersurface  (6.13)  need  not  be  flat*  This 
cc 

raises  the  question,  what  curves  on  the  hodograph  have  flat  prototypes? 

Since  a^/ar^r^  » 0,  (6.13)  will  be  flat  if  and  only  if  [3] 

(6.16)  det  jj  SXj/rp  aXj./ar2,  Bx^at,  a2Xj./art  8 1 1|  - 0 

The  product  of  (6.16)  and  det  | Vjc>  ||  yields 

(6.17)  rank  of  ||  b^  dp^/dt  ||^1 

where  b*^  ~ k Qp.  Now  let  np  fyp/dt  «*  0,  n^n^  * 1.  Then  for  some  Ih,  C^, 
and  Di,  b^  - B^n^  + ^(n^^/dt  + n^p/dt  + Di(dpp/dt)(d{x/dt ).  By 

(6.17)  the  rank  of  |J  (L  np+  dp,p/dt  ||==1.  Hence  •»  CE^,  « HEh  for  some 
C,  D,  and  non -null  E^.  Now 

(6.18)  b^  = B^n^  + Et  j^C(npdp,Y/dt  + n^  dpp/dt)  + D(fyp/dt)  (dp^/dt)] 

Suppose  and  Ei  are  linearly  independent.  Let  Fia  be  two  linearly 
independent  solutions  of  = 0.  Then  Fiebap  ■ (B^F^Jn^p,  (6.6),  and 

(6.7)  imply  (6.10),  i.e.  * ^(t)  is  a characteristic. 

Next  show  that  the  prototypes  of  both  families  of  characteristics  are 
flat.  Let  = p,a€(t)  define  one  characteristic  from  each  family  through  a 

given  point  P of  the  hodograph.  At  P,  by  (6,18) 

2 t2(e=P  ' V“o,?  > ’ (V+  V (u°.!>*'V  J * 

(-  if*1  [(dpa+1/dt)(dp2p+1/dt)  + (dp1p+1/dt)(dp2+1/dt)  ] f 
for  some  f / 0.  Then  by  (6.6)  and  (6.7) 

(6.19)  bflp  (dpo+1M)  (dpp/dt)  » 0 

^**1 

where  £ is  not  summed.  Since  these  have  non-trivial  solutions  dp.Q  /dt, 

(6.19)  implies  (6.17)  for  M-a  * V-* » 


\ 
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Now  investigate  the  relation  between  tangents  to  characteristics  and 
unit  normals  to  prototypes  x5  ■ (r,  t)  of  characteristics  p ■ uf  (t). 

By  (6.13)  » ri^Bx^/ar  = 0.  For  some  5 C 

(6.20)  n £ - A* 

Since  n^dx^/dt  • 0,  then  by  (6.18) 

(6.21)  a‘  [(rT  - AY)bJp  * k>ap]  - 0 

If 

(6.22)  det|cYbJp  ♦ ||  / 0 

where  c ■ r - A , then  by  (6.19)  the  general  solution  of  (6,21)  is  a A*« 

Y y y a 

^ ,*j 

gdp,a  /dt  for  some  g.  By  (6.20) 

(6.23)  n^=  g (dUj./dyLa ) ( dpa ) ( dp^’V dt ) 

i.e.  the  tangent  to  na  ■ |^+1(t)  is  normal  to  the  prototype  of  p.^**  . 

These  considerations  and  elementary  calculation  yield 

Theorem  6.3:  If  for  a double  wave  the  rank  of  llbaf>'  bap<  Vp  II-  2- iSfE 

(1)  The  characteristics  are  the  only  curves  on  the  hodograph  with  flat 
prototypes. 


(2)  The  tangent  at  any  point  P of  a characteristic  is  parallel  to  the  normal, 
at  points  on  the  prototype  of  P,  'to  the  prototTOe’liypersiirface  "oT" the"' other  """"" 
cHarac'teri'sti'c  "through  P.  “ 


(3)  The  prototype  of  a characteristic  is  a (sound  wave)  characteristic 
hyper  surface.  “ ” 


For  the  omitted  cases,  first  suppose  (6,22)  is  false  for  all  c . 
€ V 

Then  V " I^BaBp  and  k ■ KB^Bp,  for  some  Bni  K^,  and  K.  Now  the 

Riemann  tensor  , of  u^. 


(6.2U) 

[3], 


>0  P 


°c  ,8 

is 


R 


'apy6 


b«Y  Vs 


0a6  °Py 


where tf  is  summed.  But  now  R 


chosen  so  that 


’«Py6 


0,  so  the  parameters  p,n  can  be 


(6.25)  gap  *■  8ap 
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which  implies 

(6.26)  (aVV>(s2V3vV  " 0 

so  by  (lull),  (6.2),  (6.25>)j  and  (6.26)  92  K Ba  by 

(6.?)  9^c/3na9jip  - KBaBp.  Hence  the  rank  of  ||  x^Uj./tyi^p  + 

82k^a9|ip  j|=l,  contradicting  (It. 9). 

Next  suppose  the  rank  of  ||  b^0,  b2^,  k a ^ ||<£  2.  Then 


(6.27) 


ap 


E.  a m 
x a,3 


far  some  and  a^.  If  EQEa  * 0,  then 
(6-28)  * vi  £ - 0 

which  implies  ® 0.  Choose  parameters  such  that  (6*25)  and  (6.26) 

hold.  Then  9 Uj./aua  Qjj.p  ■ 0,  and 

(6.29)  Uj.  = aIat,.a  + bj 

where  and  bj.  are  constants  with  a^  a^  “ 6cp»  Since  vI4  may  be 

chosen  to  have  constant  components,  the  hodograph  is  the  intersection 
of  the  two  hyperplanes  v^Uj.  - b^)  = 0.  This  implies  that  the  velocity 

potential  must  be  of  the  form 

(6.30)  4 - bj-Xj.  + F(y1,  ifJ2) 

where \pa  » a^  x^  and  F is  arbitrary.  If  both  a^  **  0,  the  flow  i s a 

steady  plane  flow  in  planes  normal  to  the  lines  a^a(x^  - x^).“  0.  Now 

assume  a is  non-null.  Since  the  rank  of  a_  is  two,  rotate  space  axes 

to  obtain  the  same  form  for  $ with  - t cos  a + x^  sin  o » x2  sin  p + 

(x^cos  a - t sin  a)  cos  p for  constants  a and  p,  and  suitably  redefined 

bj..  Since  ?■<$/ u^  ■ by  these  are  swept  x^x2  - plane  flows.  If 

sin  a sin  p / 0 the  basic  plane  flow  is  any  plane  flow  steady  relative 
to  axes  translated  with  uniform  velocity  x^,  x}>,  0 defined  by  cos  a + 

x^  sin  a ■ 0,  (x^  cos  a - sin  a)cos  p + x£  sin  p - 0.  Now  |ia  * 9F/kya. 
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If  sin  or  sin  p « 0 ,"hy  a rotation  of  x-x2  - axes  at  most,  $ can  be  reduced 
to  the  form  ) - h^Xj.  + F(x1,  t),  which  corresponds  to  any  swept  unsteady 

one-dimensional  flow. 

Finally,  suppose  0.  Then  c€  b^  * 0 for  non-trivial  solutions 

of  ceEe»  0.  Since  uQ  ■ constant  leads  to  steady  flows,  assume  ^ ■ uQ, 

H2  * u3«  Then  c€v  32u  Ai,aa*ip  " °*  If  both  ce\£m  °>  then  \em 

for  some  r^.  Now  vQ£  - E^  ra9  UC/9P]_  and-  v^»  - E£ra9ua/9tJL2*  ^k051  ®6€* 

VI6VT£  * E6E£  [(raauaSuaAii)2  + rara  +-  (radna/dfi2)?J impossible.  Hence 
assume 

(6.31)  a2u2/^na9|i,p  * cs^/a^sup 

If  c is  a constant,  then  by  (6.31)  the  hodograph  lies  entirely  in  a hyper- 
* 0 for  some  constants  a^  and 

a rotation  of  space  axes  this  yields  dtf/dXj  - Adtf/dt  • B for  some  constants 

A and  B,  with  the  general  solution  ff  » Bx^  + F(Xp  Xg,  Ax^  + t).  If  A » 0 

this  is  a swept  unsteady  plane  flow,  which  still  has  to  be  chosen  to  have  a 
two-dimensional  hodograph.  If  A / 0,  the  flow  will  be  steady  relative  to 
axes  translated  with  uniform  velocity  0,  0,  l/A.-  If  c is  not  constant,  the 
integrability  conditions  of  (6.31),  (d  c/Bfi^ ) (9  u^/a^ap.^)  3 

(Bc/app)(a2u1/aMaapY)  yield  a2HrAittanp  “ a^Bp  for  some  Ba  - 9c/9|ia  and 

some  a , whence  b^;  = B„  BQ  for  some  , If  B_  B„  / 0,  then  by 

Y*  ctp  tap  o ct  & 

(6.27)  k R = KB  Br  for  some  K,  which  yields  a previously  obtained  contra- 

i®r  Op  ^ 

dicition.  If  B B » 0,  then  b . = 0,  a case  discussed  in  the 

^ a a ap 

preceding  paragraph. 

The  foregoing  discussion  and  an  application  of  (6.17)  yield 
Theorem  6.li:  If  for  a double  wave  the  rank  of  ||  b^ , b«p»  then 

(1)  The  flow  is  a uniformly  translated  swept  steady  plane  flow  or  a swept 

unsteady  one-dBmensional  flow.  ~ 

(2)  Any  curve  on  the  hodograph  has  a flat  prototype. 

(3)  Conclusions  (2)  and  (3)  of  Theorem  6.3  apply  to  these  flows. 


with  a.,  a.. 


/ 0.  After 


plane  a^Uj-  - b^) 
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7.  Triple  waves.  When  n - 3,  let  g.  • be  the  covariant  metric  tensor  of 
the  hyper  surface"  Uj.  ■ Uj(p.)  based  on  J the  euclidean  metric  ds1-  ■ du^du^., 

and  let  b.  . be  its  second  fundamental  tensor  [_3J  By  definition 

(7.1)  gi;j  - (SuI/9ni)(3u;i;/a^;j) 

(7.2)  b^  - dpi 

where  is  the  unit  normal  defined  to  within  sign  by  (lull).  Also 

(7.3)  v^a^/ap.^  q3q/3{Xi 

where  q2  - Write  aq/9^i  * q where  the  subscript^  denotes  the 

covariant  derivative  with  respect  to  p.^  and  based  on  g^. 


By  (U.l5),  (U.16),  and  (7.1) 

(7.U)  *io  - - Ok/a|Jti)gij(aur^ljij) 


x^here  gM  is  the  inverse  of  g^.  Since  the  Ghristoffel  symbols  of  the 
first  kind,  based  on  g.  are  [ij,  k]  « (8^/3^  ) (3^ ^ ) , the  second 

covariant  derivative  of  k becomes 

(7.5)  k a2k/3p.8p. . - (ak/apth)ghk(au1/ap,k)(a;?u:!:/a,iiaix.) 

^ J*J  v ***  ^ 


Now  (U.17)  becomes  (k^.  ♦ rb._  j )A jk  » - 6^,  where  r-rr  Since  (U.18)  is 

homogeneous  in  , only  the 

(7.6)  adjoint  of  (k^  + rb^)  - r”1"3!^ 


is  required,  where 

®ij  * R(k,i3' 


k,u)' 


Eij 


R(k,ij»  bij^ 


2 *i  j 


R(bij*  bij^ 


(7.7)  R(k,ij»  bij^  “ k,i+lj+lbi+2  j+2  “k,  i+ljt2bi+2j+l 

+ k,i+2j+2bi+l+lj+l  " k,i+2j+lbi+lj+2 
Since  all  four  of  r and  ^ are  independent,  (U.18 ),  (7.6),  and  (7.7)  yield 

(7.8)  [a2(glj  - u0  >1%j3)  - (»0)i*  « !><».,,  * M j)] B?,  - 0 

Since  is  a relative  contravariant  tensor  of  weight  , the 
system  (7,8)  may  be  considered  to  be  tensor  equations.  Furthermore,  note 
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that 


(?.9)  %j  - bi+lj+lbi+2j+2”bi+lj+2bi+2j+l  “Ri+li+2j+l;j+2 

are  the  components  of  the  Riemann  tensor  of  the  hodograph  hypersurface. 

(7.8)  is  a system  of  three  second  order  partial  differential  equations 
for  the  five  functions  u^.  and  k.  If  k / (■)  0,  to  determine  u^  (p. ) and  k (p) 

requires  two  (three)  additional  equations,  which  may  be  obtained  by  assigning 
a special  form  to  the  coefficient  tensor  of  (7,8).  The  resulting  systems  are 
classified  according  to  the  nature  of  the  characteristic  surfaces  of  their 
internal  hyper  surfaces.  In  the  present  context  a characteristic  is  a surface 
on  which  the  coordinate  functions,  their  first  partial  derivatives,  and  hence 
the  metric  tensor,  are  continuous,  while  the  components  of  the  second  funda- 
mental tensor  may  have  discontinuities  there. 

One  might  hope,  by  analogy  with  Theorem  6.1,  that  the  characteristics  of 
the  hodographs  of  triple  waves  would  be  two-dimensional  hodographs.  Although 
it  has  not  been  possible  to  prove  this,  it  will  be  shown  that  two  of  the 
equations  far  these  characteristics  are  consequences  of  the  equations  (6.6) 
and  (6.7)  for  two-dimensional  hodographs. 

To  find  the  characteristic  conditions  corresponding  to  (7.8),  con si der  a 
triple  of  functions  NL  (tp  t?)  of  class  with  Jacobian  matrix  ||  afh/dt^H  of 

rank  two.  Assume  that  on  the  surface 

(7.10)  p±  * Mi(t1,  t2) 


the  values  of 

(7.11)  \i±  (M(t»  - U(t),  k(M(t))  - K(t) 

(7.12)  dUjAu  * PIi(t),  dk/fy^  » Qi(t) 

are  known.  Also  on  this  surface  let  « 


and  k,.  . ■ K, 


°i y bij  ' "ij  ' 'xj 


be  the  values  of  the  three  dimensional  metric  tensor,  second  fundamental  tensor, 
and  covariant  derivative  of  k,  which  will  now  be  calculated.  Define  non -trivial 
Mj,  and  Nj  by 

(7.13)  Nj  rTi-  0,  Nj  Nj  » 1 

(7.1*0  Nj.y  BUj.Ma  - (NIr  Pn)  m±/Qta  = 0,a,y  - 1,  2. 

0 • 

- 0 


NI1NI2 


(7.1?)  Lj  SMjMa 
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By  (7*13)  and  (7.1U) 


(7.16)  Nj  »■  G€  Mi€  € ® l,  2 

for  sane  (t).  Since 

P '' 

(7.17)  S * Pu  z\/?\  #t8 

then  by  (7.2 ) 

(7.18)  Bi;j(9Mi/3ta)  (aMjMp)  “ NjB5  Uj/St^  3tp 

Also  (7.13)  implies 

(T.M)  Bi;j  ™3/ato  = - r^.  a^Ma 

Since  an/cita  and  are  independent 

(7.20)  B..  = Hap  (aMi/ata)  (aM./tp)  + Ha  (LjaMi/ata  + i±  8Mj/atc>  bl^ 


for  some  H Q H . and  B.  Let 
ap,  <x 

(7.2i)  hap  - (a-yay  (ai-i  /atp), 

(7-22)  *>"  - 32  V»a  9tp 

Then  by  (7.20),  (7.18),  and  (7.19) 

(7.23)  H^p  - 06  b«  h“V  hf* 

(7.21.)  fit,I.1V,Wh|Sr  “V*.  * hi  * 0 


Since  these  determine  Hflp,  and  Hfl,  everything  in  the  right  hand  member  of 
(7.20)  except  B can  be  found  in  terms  of  surface  data  on  (7.10). 


Next,  consider  K,^.  first  note  that 
(7.25)  3PIi/ato  - (a^y&^a^)  (3M./ata) 
By  (7.5)  and  (7.25) 


hV" 

(7.26)  K.ij  S«,Ma  ■ - V 
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Now  let  Kvj  be  the  second  covariant  derivative  of  K(t)  with  respect  to 
ta,  tp  based  on  g*‘  p - (SUj/b^)  (BUj/atp). 

By  (6.5) 

(7.27)  K\p  - a2K/ata  9tp  - OK^tY)g^6( aUj/stg )b\/dtadt? 

By  (7.5),  (7.25),  and  (7.27) 

(7.28)  K,  ±.  (91^/9^)  (3M./3tp)  - K*j  ap  - 

- <¥=%>  - */*y  g*  v8aUlM6)  3 Vv*p 

Since  the  vector  in  parentheses  in  the  right  hand  member  is  normal  to  ai^/kt^, 
then  for  some  Sc  (t)  QkGkhPih  - BK/at^g*  ^Bl^/Stg  - S£  Nj.£ 

(7.29)  K.y  OKj/St^)  - K*!  ^ - S^b*,, 

As  in  the  case  of  B.  . , expand 

J-J 

(7.30)  K,i j «Map(aM/at0)  (3mMp)  + Ma(LjaMi/ata  + i^/av  GLiLj 
for  some  v*  Ma,  and  C.  Then  from 

M«pV  V “ x>13  3Ml/3VM3M6 
h\  haY  * Mcj)  V - *«« 

and  with  the  aid  of  (7.25)  to  (7,29)  one  can  find  M&p  and  Mc  in  terms  of 
surface  data. 

Only  B and  0 remain  to  be  found,  if  possible,  from  (U.18)  which  may 
be  rewritten  as 


Fij  adjoint  (K,^  + rB..)  = 0 


for  all  r,  where  B.  . and  K,.  . arc  to  be  evaluated  by  means  of. (7.20)  and. (7. 30), 

J -*•  J 

where 


(7.31) 


Pij 


(Gij  " Po»  iPo,j> 


" <Po,i  + Vhi>  <P 


o,j  + UkPkj) 


V 
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The  fundamental  equation  may  also  be  rewritten  in  the  form 


(7.32)  fJ.  adjoint  (K1,  ±i  + B^)  - 0 
where  F^  = F±.  Jk;j,  K1,^  + r B1^ 


(K’ij  + 


*V 


JhiJkj» 


and  ||  ||*  ||  9Mj/&t2,  jj  . Now  3 and  C appear  in  the  adjoint 

‘2 

in  (7.32)  only  in  the  combination  (C  + r B)  (L^)  . Then  B and  C will  fail 
to  be  uniquely  determined  if  the  coefficient  of  this  product  vanishes,  i.e,  if 


*"«  (,(-!)“  * P 1>!  „ * ! p * ! * (P  C£  - S€  )»f*  If  *l] 


where 

ft 


*2-/  * * *0 

a p V8ct  p o j a o * p 


-<”V  * «*,  <v , ♦ tV,  p) 


Since  r is  arbitrary,  this  yields 
(7.33)  0€Pip  urfb^((1.0 

(7.31.)  p;  p (-1)“  *%*  „ * j.  p * !-Se  »a  * 1>  - 0 


Owing  to  the  presence  of  C t Sg  (7.3:0  and  ( 7.3U)  are  implied  by  but 
not  equivalent  to  (6.6)  and  (6.7).  The  three  consistency  conditions 
obtained  from  (7.33)  by  requiring  the  coefficients-  of  1,  r,  r2  to  vanish 
and  formally  setting  B * C » 0 may  conceivably  contain  additional 
information  which  can  be  combined  with  (7.33)  and  (7.3k)  to  obtain  (6.6) 
and  (6.7) i but  the  author  has  been  unable  to  reach  a conclusion  on  this 
matter. 

It  should  be  remarked  that  the  form  of  (7.33)  can  be  simplified 
considerably.  By  (7.2),  (7.13),  (7.16),  and  (7.22)  C£b*|  * N^Uj-ykt^tp 

■ vw  where  = 9M^/0ta.  Then  (7.33)  becomes 

Wij(-  vywvi  • 0 

Since  the  terms  with  r * i and  s * j vanish,  this  may  be  rewritten  as 
(7.35)  (Bi/i+1j+1  - Bi+i/ij+i  “ Bij+iFi+ij  * Bi+lj+lFij>Li+2Lj+2  * 0 


\ 
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where  solutions  of  (7.15)  way  be  defined  as 
(7.36)  l1+2  - 

In  a given  triple  watve  B.  .(M)  and  F.  ,(M)  are  known  functions  of  M.  • If 

X j 1 J 1 

the  form  (7.35)  is  not  definite,  then  as  one  would  expect,  (7.35)  with 
defined  by  (7.36)  is  a first  order  partial  differential  equation  for  the 
Mi^l*  ^2^*  ^ k a °>  as  an  nex^  (7.34)  is  satisfied 

trivially.  If  k i 0,  one  must  also  verify  that  (7.34)  is  satisfied 
oy  the  selected  solution  of  (7.35). 

8.  Examples.  Homogeneous  Flows: 


To  construct  a simple  wave  k (p ) may  be  chosen  arbitrarily.  To  construct 
double  br  triple  waves  it  must  be  chosen  to  satisfy  (6.7)  or  the  first  two 
equations  of  (7.8),  In  all  cases  k » 0 is  a possible  solution,  and  then  by 
(4.7)  the  corresponding  p are  homogeneous  of  degree  zero  in  the  x_.  This 

suggests  considering  flows,  not  necessarily  isentropic,  in  which  u. , p,  p, 

and  s are  homogeneous  of  degree  zero  in  the  x^.  If  all  of  these  functions  are 

independent  of  t,  a steady  conical  flow  field  is  obtained.  If  some  of  these 
functions  actually  depend  on  t,  let  (0.1)  X.  » x. /t.  Then  u.  , p,  p,  and  s may 

be  taken  to  be  functions  of  X. , and  the  geometry  of  the  flow  field  simply  ex- 
pands at  a uniform  rate  with  the  passage  of  time.  Let  * d/dXp  B/BXg, 

d/dly  Then  (2.1  v)  to  (2.4  v)  become 

(8.2)  (u-Y).  V*u«  -p"1  V*P 

(8.3)  (u  - Y).  V*  p + pV*n  « 0 

(8.4)  (u-Y^Vs-0 

3 

In  a homogeneous  flow  it  is  natural  to  seek  characteristic  hyper surfaces, 
defined  by  f(X)  - constant.  These  are  conical  hyper surfaces  with  vertices  at 
the  origin  of  space-time.  Now  the  char acteri stic  conditions  take  the  form 

(8.*5)  (u-xj.V*  f-0 

0.6)  [(u-Y).V*  f]2«a2(V*f)2 

Note  that  if  u is  too  close  to  X,  for  example  u * Y,  there  are  no  solutions  of 
(8.6),  i.e.  one  cannot  construct  a conical  characteristic  hypersurface  of  the 
sound-wave  type. 


\ 
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Three  classical  one  dimensional  examples  of  degenerate  unsteady  flows 
come  immediately  to  mind.  First,  the  early  stages  flow  in  a shock  tube. 

Second,  the  very  closely  associated  flows  produced  when  a piston  impulsively 
starts  to  advance  into  or  recede  from  stagnant  gas  at  constant  speed.  Third, 
the  early  stages  of  Lagrange’s  problem  in  interior  ballistics,  in  which  a 
mass  of  gas  initially  in  a uniform  state  of  rest  propels  a projectile  down 
a tube  without  friction.  These  suggest  the  following  generalizations  to 
homogeneous  unsteady  flows, 

1.  Flows  about  cones  in  a shock  tube  of  infinite  cross-section. 

Consider  flows  produced  when  an  infinite  plane  shock  moving  with  uniform 
velocity  into  stagnant  gas  and  followed  by  an  infinite  region  of  uniform  flow 
encounters  a conical  obstacle.  Such  flow  include: 

A.  Development  of  Prandtl-Meyer  flow  around  a corner  (Fig.  1A  6,  7,  8). 
This  would  also  arise  iii  the  early  stages  of  the  flow  out  of  a plane  shock 
tube  with  a flared  exit  (Tig.  IB). 

B.  Development  of  steady  conical  ilow  about  cones  not  inclined  too 
much  relative  to  the  flow  (Fig.  2) 

C.  Shock  diffraction  about  and  reflection  from  conical  obstacles  at 
large  inclination  to  the  flow  (Fig.  3)  [l, 

2.  Conical  Piston  Problems. 

These  involve  the  flow  about  a conical  object  inpulsively  started  into  uniform 
translatory  motion  advancing  into  or  receding  from  stagnant  gas. 

3.  Conical  Boundaries  in  Relative  Notion, 

An  example  of  this  type  is  an  idealized  model  of  an  emerging  projectile.  Con- 
sider a plane  block  sliding  along  a fixed  wall  with  a corner,  and  eventually 
separating  from  the  corner.  Assume  that  behind  the  block  there  is  a region  of 
uniform  flow  at  high  pressure,  while  between  the  block  and  the  wall  is  stagnant 
gas  at  low  pressure.  Consider  the  homogeneous  flow  field  after  the  block  sepa- 
rates from  the  wall  (Fig.  it). 
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